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Abstract 

The purpose of this note is to present a formulation of a given nonlin- 
ear ordinary differential equation into an equivalent system of linear 
ordinary differential equations. It is evident that the easiness of a such 
procedure would be able to open a new way in order to calculate or 
approximate the solution of an ordinary differential equation. Some 
examples are presented. 
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1 Introduction 

Let us consider an initial value problem (IVP) of first order like 



ft = 



y{to 



(1) 



or 



^ dy 
dt 



2/0, 



v(y) + 



(2) 



where the function ip{y) is supposed to be sufficiently smooth on an open 
interval (to — + 5) for some (5 > 0. 

Usually, problems like ([I]) or ([2]) are solved by means of known special 
methods according to the particular type of function for example, separa- 
tion of variables, homogeneous equations, linear equations, Bernoulli equa- 
tions, etc. However, the difficulty to find out an exact solution is related 
to some particular facts associated to the particular method of resolution. 
For instance, when the method of separation of variables is attempted to be 

applied the existence of a general primitive function of — is required; but 

the fact of computing such a primitive is not a trivial challenge. 

The goal in this work consists to transform problems like ([1]) or ([2]) into 
a following system of infinite linear ordinary differential equations 

X' = AX(t), or X' = AX(t) + b(t) (3) 



such that X(t) 



xi{t) 

X2{t) 



V 



in an infinite dimensional function vector, 





aoi 


ao2 


O03 




an 


ai2 


Ol3 


0.20 


021 


0.22 


023 


«30 


031 


032 


033 



sional matrix and 



is a real constant entries infinite dimen- 
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/ 9{t) \ 



h{t) 









is an other infinite dimensional function vector. 



v 



/ 



The organization of this work is as follows. The aim of Section 2 con- 
sists to display some particular transforms through some simple examples. 
Section 3 is devoted to expose the procedure to transform an initial value 
problem into a system of linear ordinary differential equations. In this sec- 
tion the equivalence between the solutions of the IVP and the system is 
shown. In Section 4 some suggestions about possible generalizations are 
presented. 

2 Some examples 

Before developing the details of the procedure, two examples are presented 
in order to show how this technique works. 

Example 2.1. Let us present an easy example because it is an initial value 
problem where the ordinary differential equation is linear. That is, 



for some real number yo- 

• It is straightforward to compute the exact solution because the ODE 
is linear. Then the general solution is 



for an arbitrary real constant K. Using the initial condition the con- 
stant \s K = yQ — 1. So, the solution of the IVP is 





y{t) = KJ + t + l 



y = {yQ-l)J + t + l. 



• Let us present the transformation putting 

XQ{t) = y{t) 



Oscar A. Barraza 



then 



So, 



^oi^) = y'{t) = y{t) - 1 = xi{t) :. x'q = xi. 



x'i{t) = {y{t) - t)' = y'{t) - 1 = y{t) -t-l:= X2{t) :. x'^ = X2. 
One more step 

4(0 = iy{t) -t-iy = y'{t) - 1 = y{t) -t-l:= X3{t) :. x'2 = X3. 
Then, in general for all integer j > it is deduced that 

Xj — -^j + i ■ 
In this way the matrix system is 



X'(t) 





1 





• 


" ^ 




/ 


-t \ 








1 


• 




















1 • 




X(t) + 
















• 









V : 












V 





= A = b{t) 

It is recalled that the exponential matrix of a matrix B belonging to 
a Banach algebra A is defined by 



I + B + iB2 + ^B3 + ... 



+00 



(5) 



lim S„ 

hoo 



n-1 



i=i 

where I is the identity matrix and 5"^ = I + y^ ^B-'' is the nth-partial 

j=i ■'' 

sum. This series is absolutely convergent in the Banach algebra A 
because on one side 



IB^^II < llBf 



for all integer A; > 1, and then 
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showing that the exponential matrix is bounded in the norm of 
the Banach algebra A for all matrix B A. On the other side the 
sequence of the partial sums {Snj^ is Cauchy convergent. For the 
details the reader is referred, for instance, to ([5|). 

Taking into account the variation of constants formula, the solution of 
the previous linear system is 

X(t) = e^-*.X(0) + [\^'^'-'\b{s) ds 
Jo 

and, as a consequence, the solution of the given IPV is then 



y(t) = xoit) = (1st. row of X(t)) 

/ CO \ 

Cl 

(1st. row of e^*). C2 

C3 



+ / (1st. row of e'^^*-")] 
Jo 



1 



1 







V : / 



ds 



Co + cit + —C2t^ + —cst-^ + ■ ■ ■ + / {-s)ds 



2! 
1 



3! 

+2,1 ^„+3 



= Co + Cit + —C2t^ + —C3t-' + ...-- 

1 , 1 o t 
= Co + Cit + -C2t^ + ^Csf^ + 7 

This implies that 

xo{0) = Co = yo 
and that from the equations of the system 



so, the solution is written as 
1 



for J > 1 and j ^ 2, 



C2 - 1, for j 



y{t) =yo + yot + 



2! 



yo 



(yo-l).e* + t + l. 



\jt^ + ^_{yo-i)t^ + --- + ^{yo-i)t^ + ... 



• The reader can observe that through both procedures the same solu- 
tion is reached for all real t. 
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Example 2.2. 

\y = M = y (6) 

I y(0) = yo, 

• First of all, the general solution is presented. This solution can be 
easily computed by means of the separation of variables method. The 
initial value yo must be different to zero in order to guarantee the 
nontrivial solution. Then, under the assumption yo 7^ 0, the exact 
solution of problem ^ defined for t 7^ ^ is given by 

^-yo{t-to) i-yo-t 

since to = 0. It is mandatory to keep the domain 

- — ,+00) if yo > 
yo J 

or 

00, - — r ) if yo < 0. 

\yo\J 

For yo = the solution obtained is the trivial one y{t) = 0, for all t. 

Now, we show how to transform problem ([6]) into a system of infinite 
linear ordinary equations. 

Let us define 

xo{t) = y{t) 

with xo(0) = yo. Then, 

:=x'Q{t) = y'{t) = y\t), 

so 

x[{t)=2.y{t).y'{t) = 2y\t). 
The next variable X2 is defined by 

X2{t)=y'^{t), 

and then 

x'2{t) = 3.y\t).y'{t) = 3y\t). 
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Continuing with the precedent scheme let us define 
and then 

x'^{t) = W{t).y'{t) = 4y5(t), 

and so on. 

The general expression in the present example is 
for all integer j > 0. 

Thus the set of the infinite linear differential equations can be sum- 
marized by 



=j.Xj{t), 



for all integer j > 1, 



with the initial condition xo{to) = yo- So, this system is expressed in 
its matrix form as 



X'(t) 







1 


2.X2{t) 







3.X3(t) 







4.X4(t) 











V : 



/ 1 



xiit) 



\ 



J 



where 



and 



A.X(f) 



X(t) 



xi{t) 

X2{t) 

xsit) 



\ 



I 



/ 1 

2 

3 

4 

\ \ \ \ \ \ 



(8) 
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The general solution of system ([8]) has the form 

/ CO \ 

Cl 



where 



/ CO \ 

Cl 
C2 



C2 
C3 



is the vector of arbitrary constants. 



(9) 



Finally, taking into account ([9]), ([5]) and each one of the powers of the 
matrix A, that is 



/ 2 

6 

12 

20 



V : : : : 

/ 6 



3 







24 








60 




\ 



and so on, the general solution of problem 



is given by 



y{t) = xo{t) = (1st. row of e^*^ 



/ CO \ 

Cl 
C2 
C3 

V : / 



1st. row of |^I + ^-T|(Atyj 



/ CO \ 

Cl 
C2 
C3 

V : J 
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+00 ^ 

1st. row of I + ^ - (1st. row of (A)^') 



/ CO \ 

Cl 
C2 

V : J 

Co + cit + ^.2!c2t2 + ^.Slcgt^ + • • • + j^.jlcjt^ + ■■■ = 

: Co + Cit + 02*2 + Cst^ + ---+CjP + ... 



Since the initial condition satisfies yo = y{0) = xq{0) = cq and for 
all integer j > it is clear that Xj{t) = (y(t))-'"^^, it is deduced that 
Cj = Xj{0) = {y{0)y^^ = qj"*"^ for all j > 1. As a consequence the 
solution calculated above takes the form 



y{t) = Co + cit + c|]t2 + c^tS + . . . + (J+^iJ + . . . = 

= Co [1 + cot + (cot)2 + {Cotf + ■■■ + {cQty + ...]= ^^Q-^ 



1 yo 

Co- 



1-co.t l-yo-t' 
defined under the assumption |yo-*| < 1- Note that in the open interval 
- — — 7 ) solutions (17]) and (llOp agree each to other. According 

|yol |yol/ 

to the sign of yo; the domain of this solution could be extended to +00 
or to —00. 

For the special case where y^ = 0, the solution by the series expansion 
is the zero solution as before. 



Example 2.3. 




(U) 

yo, 



• As in the previous example the solution of problem ([6|) can be rapidly 
computed by means of the separation of variables method to obtain 
for all initial value yo that 



y(0 = -ln(e"J'"-t), 



(12) 
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Its domain, for yg given, is defined by the condition e — t > wliich 
is equivalent to t < e~^", that is e^".i < 1. 

• Let us pass to the transformation of problem (jlip into a system of 
infinite linear ordinary equations. 

Let us define 

xoit) = y{t) 

with xo(0) = Uq. So 

XI (t) :=x'oit) = y'{t) = eyit), 

and 

Next the variable X2 is defined by 
and then 

x'2{t) = 2.e^y^'\y'{t) = 2e^y^'l 
Let us continue defining 

thus 

and so on. 

The general expression in this third example is 

xj{t) = e^-y^'\ 

for all integer j > 1- 

Therefore the set of the infinite linear differential equations can be 
briefly expressed by 

x'j_i{t) = {j — l).Xj{t), for all integer j > 1, 
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with the initial condition xo(io) 
expressed in its matrix form as 



Uq. Consequently, this system is 



( Mt) \ 

X2{t) 

2. X3{t) 

3. X4^{t) 



/ 1 

10 

2 

3 



V 



... \ 










Xi{t) 





















A.X(i) 



where 



X(t) 



xi{t) 

X2it) 
X3{t) 



V 



/ 



and 



/ 1 

10 

2 

3 



V ■■ 



I 



The general solution of system (jlSp has the form 



X(t) 



.At 



/ CO \ 

Cl 
C2 
C3 

V ■■ J 



where 



/ CO \ 

Cl 
C2 
C3 

V ■■ J 



is the vector of arbitrary constants. 



(13) 



(14) 
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Finally, considering ()14p , ([5]) and each one of the powers of the matrix 
A, that is 



/ 1 

2 

6 

12 



/ 2 

6 

3 24 





V : : : : : : / 

and so on, the general solution of problem (llip is given by 



y{t) = xo{t) = (1st. row of e 



At\ 



/ CO \ 

Cl 
C3 



/ 1 

1st. row of I + X] ~(^*) 



/ CO \ 

Cl 
C2 
C3 



1st. row of I + ^ ^ (1st. row of (A)^) 



+ 00 



/ CO \ 

Cl 
C2 
C3 

V J 



Co + cit + ^.C2t2 + i.2!c3t3 + . . . + i.(j - l)!cjii + 

= Co + Cit + f t2 + + . . . ^ + ^ _ 
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Taking into account that the initial condition is yo = y(0) = xo(0) = cq 
and that Xj{t) = e-''^^*^ for all integer j > 1, it is deduced that cj = 
^ji^) ~ e^'^^^^ = e^'^" for all j > 1. Therefore the solution computed 
above takes the form 

y{t) =yo + e^H + ie2-y"t2 _^ igS-yotS + . . . + ^e^-vop ^ = 

= yo + (ef«.i) + ^{eyo.tf + ^(eS'o.t)^ + • • • + l{ey<^.ty + ■■■ = 
= yo - In (1 - eyo.t) = yo-ln [e^^o. (e-^o - t)] = 
= yo-yo-ln (e-J^" - t) = - In (e-J^" - t) , 

(15) 



provided that le^^.tj < 1. In the open interval (— e ^°,e solutions 
(jl2p and (|15p are equal and the last one (jlSp can be extended to the 
interval (-00,6-2^°), that is .t < 1 as in (fT2]) . 



3 Transforming an initial value problem into a sys- 
tem of linear ordinary differential equations 

The aim of this section is to establish the relationship between the solution 
of an initial value problem of kind ([T| or ([2|) and the solution of a system of 
infinite linear differential equations of kind 

Let us begin with a necessary condition to be satisfied by the solution 
X(t) of a system of infinite linear differential equations ([3]) 

X' = AX(i) + h{t) 



such that X(t) 



\ 



is an infinite dimensional real function vector, 



/ ^00 OOl ao3 

aio ail 0,12 ai3 

^20 022 023 

^30 031 032 O33 

V ; : ; : 



is a real constant entries infinite 
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dimensional matrix and b(t) 



/ 9it) \ 






is an other infinite dimensional 



real function vector as it was presented in the introduction. Let us suppose 



that a. 



a 



for < i < j < +00 and for j + 1 < i < +cx), and that 
-1 7^ 0, for all integer j > 0. In this way the matrix A has the explicit 



form 



/o 







am 














0-12 














023 














(234 



\ 



V: : : : ■. ■■ j 

and the equations of the system are clearly given by 

x'q = aoi xi + g{t) 
x'l = a\2 X2 
x' = 023 2:3 



(16) 
(17) 

(18) 



x' = aj j+i Xj+i for j > 



(19) 



Let us also assume that the function 99(1/) is real analytic at the value y = yo 
and the function g{t) is real analytic at the value t = to- Therefore, there 
exists a unique solution y{t) in some open interval centered at the value to 
for the IVP which includes the IVP 



dy_ 

dt 

yito) 



y = -rr 



yo, 



Extending the variation of constants formula for the finite dimensional case 
( ([2]) or ([3])), it is deduced in a straightforward way that the solution of 
the system is 



X(0 



Jto 



I ds 



(20) 
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where the infinite dimensional vector of the arbitrary real constants C co- 
incides to the vector of the initial values X(to), that is 



/ CO \ 

Cl 
C2 

V '■■ J 



X(to). 



Remark 3.1. Minimal hypotheses on smoothness of the vector function op- 
erator on a Banach space of infinite dimensional vector functions X(t), 
where T(X)(t) = A.X(t) + b(t) for the case in which A is a constant ma- 
trix, or T(X)(t) = A(t)X(t) + b(t) for the case in which A{t) is a function 
matrix, guarantee a unique local fixed point of operator T, and then a unique 
local solution of the former system of infinite ODE with the initial vector 
value C. Besides, as is shown in (f^), the exponential of an operator on a 
Banach algebra can be defined by symbolic calculus, and then it is dijferen- 
tiable satisfying the standard rule 




As a direct consequence, equality [2U\) provides the solution for the linear 
system 



To reach the goal of this article which is the resolution of problem ([2]), 
the following equality has to be imposed 



'^It is usually known as the Picard fixed point theorem for abstract Banach spaces; 
see, for instance, ([5]). Alternatively, it is called the Contraction Mapping Principle on 
complete metric spaces; see, (jj) or ([I]). 
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y{t) = XQ{t) 

= (1st row of e^- . 



/ CO \ 

Cl 

V '■■ ) 



+ ( (ist row of e^-(*-") 



( 9is) \ 



V ■■ J 



it - to)^ it - tof 
Co + aoi.Clit - to) + aoi.ai2-C2 h aoi-ai2-a23-C3 -, h 



+ 



.A:=0 



2! 



3! 



(21) 



and then by differentiating with respect to t it is easy to arrive to 



+ 00 



ip{y) + 9{t)=y'it)=x'oit) = Y, 



.A:=0 



whence it results that the Taylor series for the function composition ip{y{t)) 
about the value t = to is 



+c>o 



^{y{t)) = 



h=l 



h-l 
.k=0 



Ch- 



it - to) 



h-l 



ih-1) 



(22) 



Taking into account the equations (16), (17), (18), . . . , (19) of the linear 
system it is evident, for all real t in an open neighborhood of to, that 

x'oit) = aoi.xi(t) +g{t), 



so 



+ 00 



aoi aoi aoi f-' 

h=l 



'h-l 

W ak,k+i 
,fc=0 



Ch- 



it - to) 



h-l 



ih-l)\ ' 



and for all integer j > 2 that 



x'j_iit) = aj^ij.Xjit), 
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so 



Xj{t) 



-, +00 

fc=0 



h-1 
.fc=0 



Ch- 



it - to) 



h-j 



{h-jy. 



Mt)) . 



1 

fc=o 

And finally, by evaluating at t = to it is obtained from (I2ip that 



Co = yito) = yo, 



(23) 



and that 



1 



(24) 



t=to 



k=0 



for all j > 1, in coincidence with the construction of the coefficients of Taylor 



series (|22]) where 



.fc=0 



The previous set of infinite equalities allows to compute the constants 
Cj for all integer j > in terms of the function (p{y), the derivative of 
y{t) = XQ{t). In fact, the following statement has been proved. 

Theorem 3.2. Let us suppose an infinitely derivable function g{t) att = t^ 
and the existence of the function y{t) solution of problem ^ 

I y'(t) = f = mt)) + 9{t) 
I yito) = yo, 

for a function (p{y{t)) which is represented by its Taylor series about t = to 

it-toy-' 



Mt)) = E 



i-i 

Y\_ "■k,k+l 
.k=0 



■CJ. 



u - ly 
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on the open interval \t — tol < 6 defined for some real number 6 > 0. 

Then, the family of infinite functions given by the following relations for 
xo{t) = y{t) 

a^o(0 = "Piyit)) + gif) = aoi.xi(t), 

x[{t) = ai2.X2{t), 
x'^{t) = a23.X3{t), 

x'j{t) = ajj_|_i.Xj-|_i(t), for all integer j > 1 
defines an infinite dimensional function vector 



X(t) 



Xi{t) 
X2{t) 

V ■■ J 



which is the solution of the system of infinite linear ordinary differential 
equations 



X'(t) = A.X(i) + b(t) 



with the initial vector value 



X(to) 



/ CO \ 

Cl 

V '■■ ) 



Proof. The assumptions about functions g{t) and (p{y{t)) allows to establish 
that the (unique) solution y(t) of the initial value problem ([2]) on an open 
interval centered at the value t = tois infinitely differentiable there thanks to 
the application of the chain rule. In this way the existence of the constants 
Cj, for all integer j > 0, is guaranteed according to the family of formulas (j23p 
and (j24p preceding this theorem. Taking into account these computations, 
it is straightforward to conclude that the function vector X(t) is the solution 
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of system ([3]) satisfying the initial conditions X(to) 
the purpose of this result. 



/ CO \ 

Cl 
C2 

V : J 



as it was 



□ 



Theorem 3.3. Conversely, let C 



/ CO \ 

Cl 
C2 

V ■■ J 



be an infinite dimensional real 



vector and if 

1. the function g[t) is infinitely differentiable on t = tQ, 

2. the function ^{y) is infinitely differentiable for all real y such that 
|y — Col < /3 for some real number /3 > 0, 



3. the power series 



E 



■i-i 

.A;=0 



■CJ. 



U - 1)! 



converges on the open interval |t — to| < Oi for some real number a > 0, 
and 



4- if, besides, system (3\) 



X'{t) = A.X(t) + h{t) 



with initial vector value X(to) = C admits as solution the infinite 
dimensional function vector 



X(t) 



Xi{t) 
X2{t) 

V : J 



20 



Oscar A. Barraza 



then, its first component y{t) = Xo{t) is the solution of the initial value 
problem ^ on an appropriated open interval \t — tQ\ < 6 for some real 
number 6 > 0, where yo = cq. 

Proof. In order to verify this theorem, let us define y{t) = XQ{t), the first 
component of the vector X(t) solution of system ([3]) 

X'(t) = A.X(i) + b(t) 



with initial vector value X(to) 



/ CO \ 

Cl 
C2 
C3 



:= C. 



V : / 

As it was mentioned above, thanks to the variation of constants method, 
the solution of this system has the form 

X(t) = e^-(*-*»).C + /* e^-(*-^).6(s) ds, 

since A is a constant matrix and then it commutes with the exponential one 
e^\ 

Therefore, the first component y{t) = XQ{t) of vector X(t), is given by the 
first row rule; that is, 

y{t) = xo{t) 

/ CO \ / 9{s) \ 



(1st row of e^-(*-*")) . 



Cl 
C2 



+ (ist row of e^-(*-")) . 








ds 



Co + aoi-Ci[t - to) + aoi.ai2.C2 -, h aoi.ai2. 023X3 \ h 



2! 



3! 



+ 



i-i 

.k=0 



Jto 



'^J T" 



Or more briefly 



+00 

y{t) =Co + Y, 



.k=0 



it - toy 

M 



+ / g{s)ds, 
to 
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which has sense for ah real t in the interval of convergence |i — to| < ct of 
the last series, for some real number a > 0. 

First, by evaluating on to it is clear that y{to) = cq = uq. 

Next, by differentiating term to term the series with respect to the vari- 
able t in the convergence interval |t — to| < and applying the fundamental 
calculus theorem it is evident that 

+00 

y'it) =E 

By hypothesis, function g{t) is infinitely differentiable on t = Iq and also 
the last power series on the open interval |t — to| < ct. Thus, the function 
y{t) is infinitely differentiable on t = as well, so by applying the chain 
rule indefinitely to the function composition ip{y{t)) on to the identities (p^ 
are rapidly recovered. As a consequence the series 

jt-toy-' 

(J-I)! 

becomes the Taylor series of the function ip{y{t)) on an appropriated interval 
|t — to| < for some real number 6 depending on to, on a, on (3 and on 
the values of g{t) on a neighborhood of to- This fact shows that y'{t) = 
ip{y{t)) + g{t) proving that y{t) is the solution of the initial value problem 
([2]), as desired. 

□ 

Remark 3.4. According to a particular situation, the domain of the defini- 
tion of the solution concerning the power series could be extended to a larger 
interval. In fact, this was observed in the two previous examples. 



.k=0 



{t - to) 



(j - 1)! 



+ 9{t). 



+00 

E 

7=1 



-i 



n 

.fc=0 



3.1 An special extension of the method 

Let us now continue with the IVP 

1 y'W = f = Hy{t),t) = ^{y{t)).f{t) + g{t) 
[ y{to) = yo, 

which its corresponding system of infinite linear differential equations is 
supposed to be 

X' = A(t).X(t) +b(t) (25) 
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such that as before X(t) 



vector, 



xiit) 



\ 



is an infinite dimensional real function 



A{t) 



( aoi/(t) 

ai2/(t) 

a23/(t) 



V ; ; : ; 



A.f{t) (26) 



is an infinite dimensional matrix with real function entries, where the matrix 





aoi 





• 


•• \ 




/ 


9{t) \ 










• 




















023 • 




has constant entries and b(t) = 
















• 









V : 








•• / 




V 





is an other infinite dimensional real function vector. Thus, the system is 
explicitly given by 



x'q = aoi.f{t)xi +g{t) 

ai2-f{t) X2 

a23-f{t) X3 



x[ 



(27) 
(28) 

(29) 



^'i = aj,j+i-f{t) Xj+i for j > 1 



(30) 



Since the matrix A(t) = A.f{t) and the integral / A(s) ds = A. / f{s) ds 

J to ^ Jto ^ 

commute, then A{t) and the exponential matrix e^^o ^^^^ '^^ = e^' ■^^o ^^^'^ 
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commute as well, and thanks to the variation of constants formula the solu- 
tion of the previous linear system can be expressed as 

X(t) = e^*>('')'='^X(^o) + /^/>W'^^b(r)dr = 

t (31) 

Jtn 



See, for instance, ([2]) or ([3]). 

Thus, under little adaptations to this case, it is not difficult to deduced 
the corresponding extensions of the previous two theorems. 

Theorem 3.5. Let us suppose an infinitely derivable function g[t) att = to, 
the existence of the function y{t) solution of problem 

1 y'{t) = ^ = Hy{t),t) = Mt))-f{t) + 9{t) (32) 
for a function if {y{t)) which is represented by its Taylor series about t = to 



^{y{t)) = E 



.k=0 



. (t - to)'-' 



on an open interval |t — to| < ^ defined for some real number (5 > 0, and an 
integrable function f[t) on the same interval. 

Then, the family of infinite functions defined by the following relations 

xo{t) = y{t) 
x'o{t) = ^{y{t)).f{t)+g{t) = aoi./(t)xi(t), 
x[{t) = ai2.f{t)x2{t), 
4(0 = a23-f{t)x3{t), 



x'At) = ajj+i./(t) a;j4.i(t), for all integer j > 1 
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defines an infinite dimensional function vector 



xi{t) 

X2{t) 

V : J 



which is the local solution of the system of infinite linear ordinary differential 
equations [25\) 

X'{t) = A{t).X{t) + h{t) 



with the initial vector value 



X(io) 



and where A(t) is given by ([26 



Theorem 3.6. Conversely, let C 



/ CO \ 

Cl 

V '■■ J 



( CO \ 

Cl 
C2 



be an infinite dimensional real 



vector and if 

1. the function g{t) is infinitely differentiahle on t = to, 

2. the function ip{y) is infinitely differentiahle for all real y such that 
|y — Col < f3 for some real number /? > 0, 



3. the power series 



E 



.k=0 



■CJ 



{t - to) 
U - 1)! 
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converges on the open interval \t — to\ < a for some real number a > 0, 
and 



4- if, besides, system ^ 



X'(t) = A.X(t) +b(t) 



with initial vector value X(to) = C admits as solution the infinite 
dimensional function vector 



X(t) 



xi{t) 

X2{t) 
X3{t) 



V 



then, its first component y{t) = XQ^t) is the solution of the initial value 
problem (0) on an appropriated open interval |t — to| < ^ for some real 
number 6 > 0, where yo = co- 
in order to illustrate this extension, let us present the following example. 
Example 3.7. 

y (33) 



^ dt 

y(0) = yo 



• The general solution is easily computed by means of the separation of 
variables method. If the initial value yo ^ the solution is a nontrivial 
one; then, the exact solution of problem ([33|) is given by 



yit) 



2yo 



2-yo-t^' 
Its domain is given by the condition 

yo.t^ + 2 



(34) 



For 2/0 = the solution obtained is the trivial one y{t) = 0, for all t, 
which can be included in the previous formula. 
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Let us pass to the transformed of problem ()33p into a system of infinite 
linear ordinary equations. 

Let us define 

xo{t) = y{t) 

with ?/o = a;o(0) = cq. Then, 



where 



So 



where 



Then 



x'Q{t)=y'{t)=y\t).t = t.xi{t), 

xi{t)=y\t). 
x'i(t) = 2.y{t).y'{t) = 2t.y\t) = 2t.X2(t), 
X2{t)=y\t). 



x'^it) = 3.y\t).y'{t) = 3t.y\t) = 3t.X3(t), 
and so on. In this way the general expression is 

x,it)=y^+\t) = xi;\t), 

and 

x'^it) = (j + l).y'{t).y'{t) = U + = {j + 

for all integer j > 0. 

Thus, the system is expressed in its matrix form as 



X'(t) 



/ t.Xlit) \ 




1 ^ 


t 








• 


" \ 


/ 


Xo{t) \ 


2t.X2{t) 










2t 





• 






xiit) 


3t.X3(t) 













3t 


• 






X2{t) 


At.Xiit) 
















At ■ 






X3{t) 


V : ) 




V : 










• / 


V 


■■ J 



A(t)= A.t 



A(t).X(t) 



(35) 
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where 



X(t) 



xiit) 

X2it) 
X3{t) 



and 



/ 1 ••• \ 

2 

3 

4 



V : : : : : 

The general solution of system ()35p has the form 



/ CO \ 

Cl 
C2 
C3 



where 



/ CO \ 

Cl 
C2 
C3 



is the vector of arbitrary constants. 



(36) 



Finally, taking into account (j36p . ([5]), the matrix / A.s ds = A. sds 

Jo Jo 

A.— and each one of the powers of the matrix A, that is 



/ 2 

6 

12 

20 

\ \ \ \ \ \ \ 
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/ 6 

24 

3 60 



V ; ; ; ; ; ; 



and so on, the general solution of problem (|3cip is given by 




Since the initial condition satisfies i/q = y(0) = xq{0) = cq and for 
all integer j > it is clear that Xj{t) = {y{t)y~^^, it is deduced that 
Cj = Xj{0) = {y{0)y~^^ = cpQ^^ for all j > 1. As a consequence the 



A reformulation of an ODE 



29 



solution calculated above takes the form 



yo 



1 2yo 

yo-— 



(37) 
2 

defined under the assumption lyo-^^l < 2, that is < -j — r, if yo 7^ 0. 

\yo\ 

rr rr ^ 

Note that in the open interval —\ -. — : < t < . - — -, solutions (|34|) 

V |yol V m\ 

and ([57|) agree each to other and solution ([57]) accepts an extension to 

rr 

all real number t ^ i^/i — r 

V lyol 

Besides, for yo = the solution (|37p is the zero series. 



4 Some conclusions 

It is impossible to deny the importance of the role of the differential equa- 
tions in all the processes which involve a dynamical component of a model 
for continuous variables. These processes are illustrated in a large vari- 
ety of domains related, for instance, to applications to Astronomy, Biology, 
Chemistry, Economy, Engineering and Physics. 

The necessity to describe the behavior of the exact or approximated solu- 
tions of those dynamical models provides an understanding of the underlying 
processes. The fact of computing the exact solutions of a differential equa- 
tion is a desirable purpose but in some cases it is a hard task to accomplish, 
in most of them it is impossible. So, finding alternative ways to solve ordi- 
nary differential equations are always good news or, at least, it could open 
some paths to obtain approximated solutions. For example, keeping in mind 
that the solution of the system 



X'(t) = A(t).X(t) + b(t) 
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with initial vector value -'^(to) 



/ ^0 \ 

Cl 
C3 



C is given by formula ([ST 



V : / 

X{t) = e^W.X(to) + r e[^W-^Ml.b(r) dr, 

J to 

where the matrix M(t) = / A(s) ds, one possibility consists in replacing 

Jto 

the exponential matrix in that formula by its n-th partial sum. In this 
way the sequence of approximations to the desired solution is defined by 



Xo(t) = X(to) = C 



and for each integer n > 1 



Xn(t) 





[M(t) - M(s)]J 



b(s) ds. 



For the special case of a constant infinite dimensional matrix A(t) = A, the 
matrix M(t) takes the form M(t) = / Ads = {t — to)-A and then the last 
sequence is reduced to 



X„(t) 




J! 




h{s) ds, 



for each integer n > 1. 

In particular, the method presented in this note which consists to trans- 
form a nonlinear differential equation into a system of infinite linear ordinary 
differential equations has the advantage to translate the difficulty to treat 
directly with nonlinearity into the problem of the treatment of the infinite 
quantity of linear ODE with constant coefficients or simple variable ones. 
The methodology of reducing the difficult of nonlinearity keeps a certain 
analogy between the methodology of reducing a nth-order ODE to a first- 
order system of n ODE because, in general, solving an isolated linear ODE 
with constant coefficients is easier than solving a nonlinear one. 



A reformulation of an ODE 



31 



However, the transformation of a nonlinear ODE to a system of infinite 
linear ODE's is only developed here for the case of first-order differential 
equations; this is the first restriction of our method. A second limitation 
is due to the hypotheses about the smoothness of functions (p and g in 
problems ([T]) and ([2]) and the existence of the series expansion detailed in 
the theorem. The third constraint of this method is the form of the derivative 
y'{t) in the first two problems. A more general situation to be studied with 
this procedure in a next work could concern an ordinary differential equation 
like y'{t) = (p{y,t), for more general function if(y,t). 
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